Long- wavelength spin- and spin-isospin correlations in nucleon matter 
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We analyse the long-wavelength response of a normal Fermi liquid using Landau theory. We 
consider contributions from intermediate states containing one additional quasiparticle-quasihole 
pair as well as those from states containing two or more additional quasiparticle-quasihole pairs. 
For the response of an operator corresponding to a conserved quantity, we show that the behavior 
of matrix elements to states with more than one additional quasiparticle-quasihole pair at low 
excitation energies ui varies as 1/uj. It is shown how rates of processes involving transitions to two 
quasiparticle-quasihole states may be calculated in terms of the collision integral in the Landau 
transport equation for quasiparticles. 
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I. INTRODUCTION 

In liquid helium 3, the prototypical Fermi liquid, in- 
teractions between atoms are to a very good approxima- 
tion invariant under independent rotations in coordinate 
space and in spin space. This fact, together with the 
translational invariance of the system, makes for an eco- 
nomical description of the properties of the system in 
terms of Landau's theory of normal Fermi liquids. The 
reason for this is that at long wavelengths, the only ex- 
citations of importance for the low-frequency response 
are those with an added quasiparticle-quasihole pair, and 
collective modes such as zero sound, which are coher- 
ent superpositions of single quasiparticle-quasihole pairs. 
An analysis of the behavior of matrix elements of the 
long-wavelength components of the density operator to 
single-pair and multipair states was given by Pines and 
Nozieres Q, and our purpose in this paper is to extend 
the analysis to more general responses. 

From investigations of the response to probes more 
general than a scalar potential it is clear that con- 

servation laws play a key role in determining the form of 
long- wavelength matrix elements of an operator. In nu- 
clear matter, interactions are not central, due to the large 
tensor contribution coming from exchange of pions and 
to other terms that give rise to spin-orbit coupling. As 
a consequence, the total spin is not a conserved quantity 
and matrix elements to states containing more than one 
quasiparticle-quasihole pair (which we shall in future re- 
fer to simply as a pair) do not vanish at long wavelengths. 
A qualitative discussion of the effects of noncentral in- 
teractions is given in Ref. 4] and Landau theory was 
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extended to take into account noncentral interactions in 
Ref. @. 

Multipair states play an important role in the con- 
text of neutrino processes in dense matter. The rates 
of a variety of weak interaction processes can be ex- 
pressed in terms of correlation functions for nucleons or 
other constituents The reason for this simplifica- 
tion is that, due to the interactions between leptons and 
hadrons being due to weak and electromagnetic effects, 
rates may be expressed as a convolution of a hadronic 
correlation function and a leptonic one. Because corre- 
lations among leptons are relatively weak, the most dif- 
ficult part of calculating weak interaction processes in 
dense matter is the evaluation of the hadronic correla- 
tion functions. A review of calculations of rates may 
be found in Ref. Q. Processes such as bremsstrahlung 
of neutrino- antineutrino pairs in nucleon-nucleon inter- 
actions and the production of neutrinos and antineutri- 
nos by the modified Urea process are both intrinsically 
two quasiparticle, two quasihole processes, and they are 
important when the corresponding single-pair process is 
kinematically forbidden. Bremsstrahlung by creation of 
a single nucleon pair is kinematically forbidden, since the 
velocity of a quasiparticle is always less than c, and the 
direct Urea process, the single nucleon pair analog of the 
modified Urea process, is allowed in dense matter only 
for certain compositions. There are many calculations of 
rate of scattering of neutrinos by nucleons that take into 
account single-pair process es M, |j| but, as Raffelt and 
collaborators have stressed [Ty| , multipair processes can 
be important quantitatively. 

The purpose of this paper is to express the two- 
quasiparticle-quasihole contribution to the response in 
terms of the collision integral in the Landau kinetic equa- 
tion for quasiparticles. Landau theory provides a clear 
separation between long- wavelength degrees of freedom, 
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II. BASIC FORMALISM 



which are treated explicitly, and short- wavelength ones ^f*~~~ P + Q/2 _ 

whose effects are included through the use of parameters x ^L^, p - q/2 _ ~ p - q/2 

that include the effects of renormalization of matrix ele- 
ments of currents and interparticle interactions. Another 
strength of the Landau theory is that it brings out clearly 
the role played by conservation laws. For simplicity, we 

shall restrict ourselves to normal systems and we shall P + 

not consider effects of nucleon pairing and superfluidity. 



+ 




p - q/2 



In Landau theory, one limits oneself to temperatures 
low compared with the Fermi temperatures of the con- 
stituents and regards the elementary excitations as being 
quasiparticles which may be characterized by an effective 
mass to*, where the index i = n, p denotes neutrons and 
protons, respectively. We shall treat the nucleons as non- 
relativistic, which is a good first approximation. 

If a system is subjected to a perturbation 



H 1 = O q c/ q e- 



Hermitean conjugate, 



(1) 



the Fourier transform of the linear response of the expec- 
tation value of the operator O q at frequency to is given 
for a system invariant under parity by [l| 

<O q >c=-X(<Z,^q, (2) 

where the linear response function is defined by 



-Ei/T 



ur- x — [w + ir})' 



(3) 



Here I denotes an initial state, j denotes an intermediate 
state and LUji = Ej — Ei , where Ej and Ei are the energies 
of the states, Z = J2je~ Ej ^ T is the partition function 
and T is the temperature. In this paper we shall only 
consider the response of the physical quantity which is 
coupled to the applied field. For systems with noncentral 
interactions, there are generally off-diagonal terms, such 
as the response of the spin in the x direction when a 
magnetic field is applied in the z direction. Since such 
effects vanish in the limit of small q, they are generally 
of less importance than the diagonal response [8( • 

A particularly useful quantity related to this is the 
dynamical structure factor, which is defined by 



jl 

= 7r (l-e-"/ T ) lmxiq,Uj) - (4) 

To make explicit calculations we shall use the language 
of diagrammatic many-body theory expressed in terms of 
fully renormalized nucleon propagator G(p,o>), where p 
is the particle momentum and lo its energy. For simplic- 
ity we shall not write out spin and isospin labels except 



FIG. 1: Diagrams for the matrix element for a transition to a 
state with one quasiparticle-quasihole pair. The dashed line 
represents the interaction with the applied external field. The 
open circle corresponds to the sum of all connected diagrams, 
and the shaded ones to all connected diagrams irreducible 
with respect to a single quasiparticle-quasihole pair. 



where necessary for clarity. We shall further write the 
propagator as the sum of a coherent part, which corre- 
sponds to a single quasiparticle or quasihole, and an inco- 
herent part, which takes into account more complicated 
intermediate states: 



G(p,Lo) = G coh (p,u) + G iac (p,Lo), 
where the coherent part has the form 



G coh (p, w ) 



u — e p + it] 



(5) 



(6) 



Here e p is the quasiparticle energy, measured with respect 
to the chemical potential for the species in question and 
z v is the wave function renormalization factor. 



III. SINGLE-PAIR PROCESSES 

To set the scene and establish notation, we consider 
transitions to states with a single quasiparticle-quasihole 
pair. The matrix element for the process is indicated 
diagrammatically in Fig. ^ by the open circle. In Fig. 
n the single particle propagator represent only the co- 
herent part. By the standard techniques of microscopic 
many-body theory this may be expressed in terms 
of diagrams irreducible with respect to states having a 
single quasiparticle-quasihole pair in the channel with 
total momentum q as shown in the figure, where the 
hatched vertices correspond to the irreducible diagams. 
Algebraically, it is given in a compact matrix notation by 



0?>=A o z(l + x/)- A = A c 



(7) 



where we suppress sums over momenta and particle 
species Here A is the irreducible vertex for interaction 
with the external field, 



3 








< 



FIG. 2: Diagrams for the matrix element for a transition to 
a state with two quasiparticle-quasihole pairs. As in Figure 
1, open circles denote the sum of all connected diagrams, and 
here the hatched circle corresponds to the sum of all diagrams 
that do not have the external field attached to one of the 
external lines. 



x P (q,w) = 



»p+q/2 - n p-q/2 
' ~ £ P+q/2 + £ p-q/2 



and 



fp P ' = z p 2y r lrr (p,p') 



(8) 



(9) 



is the Landau quasiparticle interaction, P" being the 
two-particle vertex function irreducible in the sense de- 
scribed above. We will give expressions for \° in the 
Appendix. For the long- wavelength, low-frequency re- 
sponse the only quasiparticles of importance are those 
to the Fermi surface, and we denote z PF by z. In the 
phenomenological theory, the quantity g® = zSP corre- 
sponds to an effective "charge" of a quasiparticle for the 
particular probe in question; for example in the case of 
a scalar probe, it is the total particle number carried by 
a quasiparticle and for a magnetic field it corresponds to 
the magnetic moment of the quasiparticle. For simplicity 
we have not written spin and isospin indices explicitly. 

The contribution to the imaginary part of the response 
function is therefore 

Im X sp = it Y, |O q P| 2 (n p+q/2 -n p _ q/2 )x 

p,spin 

<5(^-£p+q/2+ep-q/2) (10) 



IV. TWO-PAIR PROCESSES 

We turn now to intermediate states with two quasi- 
particles and two quasiholes. The matrix element for the 
process is illustrated in Fig. |2] Of particular interest to 
us are contributions in which the external field interacts 
directly with just one of the quasiparticles or quasiholes 
in the intermediate state. These are represented by the 



first four diagrams on the right hand side of the equation 
in Fig. |21 and they are particularly important for the re- 
sponse to operators which do not correspond to conserved 
quantities. As we shall see, they lead to divergences of 
matrix elements for low lu analogous to those found for 
bremsstrahlung. The corresponding effects in low-order 
perturbation theory have previously been drawn atten- 
tion to in Ref. 10]. The final term in Fig. [21 represents 
the remainder. We denote the quasiparticles and quasi- 
holes in the intermediate state by the indices i = 1, . . . , 4, 
as indicated. In this paper we shall consider only opera- 
tors which do not convert neutrons into protons, or vice 
versa, that is they couple only to the third component of 
the isospin. When the operator can change the isospin, 
which is the case for the response function of interest for 
the modified Urea process, the situation is more compli- 
cated, and we shall comment on this point in the Con- 
clusion. 

For small q, the contribution to the on-shell matrix 
element coming from the first term may be evaluated 
straightforwardly, and it is given by 

— ^— Af 



Qtv = _ z 2 



c pi-q 

\0 



-A- 



where 



A 



(11) 



(12) 



is the scattering matrix for two-quasiparticle scattering. 
Therefore the contribution to the matrix element from 
the first four terms is 



where 



e° = 



(13) 



(14) 



1 — v, • q/w 

The imaginary part of the two-pair contribution to the 
response function is therefore 



1234 



[nin 2 (l - 713) (1 - 714) - (1 - ni)(l - 712)71377,4] x 
S((J + ei + e 2 - e 3 - £4)6(q + Pi + P2 P3 P4)} ■ 

(15) 

where sums over quasiparticle states include sums over 
spins and isospins and the factor of 1/4 is necessary to 
avoid overcounting. Inserting the expression (|13Jl . one 
finds [13 

Wp = \Af + - - )) x 

1234 ^ ^ ' 

[7ll7l2(l — 7J 3 )(1 — 714) — (1 — 7li)(l — 712)713714] X 

S(oj + ei + e 2 - e 3 - e 4 )<S(q + pi + p 2 p 3 P4) } 

(16) 
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Let us now consider the limit u> » vpq. For external 
fields that are independent of the particle momentum, 
the quantities A are independent of the direction of the 
quasiparticle momentum. For long- wavelength probes, 
the q in the momentum conservation condition may be 
ignored to leading order, and therefore Eq. Ijl6(l becomes 

w p = £ E {M a tf(tf +tf - - ef) * 



where 



1234 



[nin 2 (l - n 3 )(l - ra 4 ) - (1 - ni)(l - n 2 )n 3 n 4 ] x 
S(w + ei + e 2 - £3 - £4)<5(pi + P2 P3 - P4)} 



(17) 



For convenience, we adopt the usual convention that the 
phase of single-particle states is chosen so that X° is real. 

It is instructive to rewrite this expression to show 
that the quantity that enters is very closely related to 
quasiparticle relaxation time. The collision term for two- 
quasiparticle scattering is given by 

UL 234 

[nin 2 {l - n 3 )(l - 714) - (1 - m)(l - n 2 )mn4\ x 
5(u> + ei + e 2 - e 3 - e 4 )<5(pi + P2 P3 P4)} , (18) 

where in this equation the distribution function is not 
necessarily the equilibrium one. On linearizing this equa- 
tion about the equilibrium distribution for quasiparticles 
"ith momentum pi, one hnds 



wi 



dn\ Sni 



(19) 



where 



— = 2tt V {|^| 2 [n 2 (l - n 3 )(l - n 4 ) + (1 - n 2 )n 3 n 4 ] : 



234 



8{uj + ei + e 2 - e 3 - e 4 )<5(pi + p 2 P3 P4)}(20) 

We now rewrite Eq. I|17|) so that the combinations of 
distribution functions is similar to that in Eq. i[2()[). By 
virtue of the relationship 



n(e) = e~ e/T [l - n(e 



(21) 



we may write 



[nin 2 (l - n 3 )(l - n 4 ) - (1 - ni)(l - n 2 )n 3 n 4 ] x 
£(cj + ei + e 2 - e 3 - e 4 ) = 
[ni - n(ei + w)][n 2 (l - n 3 )(l - n 4 ) + (1 - n 2 )n 3 n 4 ] x 

(5(w + ei + £2 - £3 - £4), 



(22) 



from which it follows that we may rewrite Eq. (|17|) as 

w p ~E(^) 2[ni ~ n(ei+w)] 



7f (ei+w) 



(23) 



2 - Ei- 4 



A? 



A 3 



A 4 ° 



Af 



1 v 1 y 234 

[n 2 (l - n 3 )(l - n 4 ) + (1 - n 2 )n 3 n 4 ] x 
5(uj + ei + e 2 - e 3 - e 4 )5(q + p 4 + p 2 - p 3 - p 4 )| 



(24) 



is the relaxation rate for the quantity X° for a sys- 
tem close to equilibrium. This gives the collision rate 
weighted by the change in A of all the quasiparticles in a 
collision. 

Depending on the properties of the operator O, the 
matrix elements to two pair states can have different 
sorts of behavior. Let us begin by considering the 
limit lj 3> VFq, where vp is a typical Fermi veloc- 
ity. The matrix element then becomes proportional to 
[A°(pi) + A°(p 2 ) - A°(p 3 ) - X°( Pi )}/uj. If the operator 
O is the particle density, X° becomes a constant at for 
u) 3> v-pq. Since the interaction conserves particle num- 
ber, it also conserves the number of quasiparticles and 
therefore in this limit , X° (pi)+A° (p 2 ) - X° (p 3 ) - A° (p 4 ) 
vanishes. It the operator is that for a component of the 
spin, a similar argument will apply if the interactions be- 
tween components arc central, that is their spin structure 
of the form of the unit operator or of the spin-exchange 
operator. However, when there are non-central interac- 
tions, e.g. spin-orbit or tensor forces, the combination of 
A's will not generally vanish. For example, two quasipar- 
ticles with spin up can scatter to states with two quasi- 
particles with spin down, in which case the matrix el- 
ement to a two-pair state will diverge as 1/lu. For the 
operator for the number density of particles, the lead- 
ing contribution to the matrix element for a transition 
to a two-pair state will have a factor v-pq/uj 2 . The ma- 
trix element will then vanish at least as rapidly as q for 
0, in accordance with the general arguments in Ref. 
lj. However, if the corresponding current is conserved 
in the long- wavelength limit, as it will be if the system is 
translationally invariant, the term proportional to q will 
also vanish, and the leading term will be or order q 2 . For 
a non-relativistic system with a single component, the 
particle current is proportional to the total momentum, 
which is conserved for a translationally invariant system. 
However, for a multicomponent system, the total particle 
current is proportional to the total momentum if all com- 
ponents have the same mass, and therefore for arbitrary 
masses of the components, the two-pair matrix element 
will be proportional to q. For a system with components 
that do not all have the same mass, it may be seen from 
this argument that matrix elements of the mass density 
operator to two-pair states will vary as q 2 for small q, 
since the total mass current is conserved. 
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V. TWO-PAIR CONTRIBUTIONS TO 
STRUCTURE FUNCTIONS 

We now evaluate the sums over momenta in Eq. (|24|l . 
For definiteness, we shall assume that quasiparticles 1 
and 3 are the same species i. Since collisions conserve 
the numbers of neutrons and of protons, this implies that 
quasiparticles 2 and 4 are the same, but not necessarily 
the same as 1 and 3. All quasiparticles participating 
in a collision at low excitation energies lie close to the 
respective Fermi surfaces, so the collision geometry may 
be specified in terms of the angle 9 between p x and p 2 , 
and (f>, the angle between the plane containing pi and p 2 
and the plane containing p 3 and p 4 . We use the result 

"GJ 



P2P3P4 

m*(m*) 2 



V £ (3) ( P i +p 2 -p 3 -p 4 ) 



m*Y f 

—7 — / 0ij sin 9d9d(j)d(j)2de2d€3d€4, (25) 
(2irf J 

where <p2 is the azimuthal angle of p2 with respect to pi 
and 



Pj 



(Pi + p] + ZPiPj cos 9) 1 / 2 ' 



(26) 



Provided the \i may be regarded as independent of the 
magnitude of the quasiparticles momentum, which is a 
good approximation for the response considered in this 
article, the integrals over energies and over momenta de- 
couple, and from Eq. (|T7Jl we find that 



t P 



1 x „ f ( m*m*) 2 



E 



spins 

J deide 2 de 3 de i S(uj + e\ + e 2 - £3 — e 4 ) x 
[nin 2 (l - n 3 )(l - n 4 ) - (1 - ni)(l - n 2 )n 3 n 4 ] j , 



where 



7y - PiPij — / ■> > 



PiPj 



(pj +P 2 j+ ZpiPj cos Q) 1 / 2 ' 



(27) 
(28) 



is symmetric: 7y = 7^;. The integral over energies is |3J, 
p. 114] 

J dt\dt2dt-ide^(yj + e± + e 2 — e 3 — e 4 ) x 



[nina(l - n. 3 )(l - n 4 ) - (1 - rti)(l - n 2 )n 3 n 4 ] = 

w[w 2 + (2^T) 2 ] 



(29) 



and therefore we find 



Imx 



tp 



(2ttT) 2 



6u> 



E 



(m*m*) 2 



spins 



(30) 



which for >> Upg reduces to, 



cj 2 + (2^T) 2 ] ^ J (m*m*) 2 
6lu ^ I (2tt) 7 



dfld4>2jij 



l^| 2 Af(Ag. + A| 



A 4j ) 



spins 



(31) 



where from Eq. (|14|l one can see that £ — > A in this 
regime. This result shows that at zero temperature, 
Imx tp is proportional to u) and independent of q. This be- 
havior should be contrasted with that of the single-pair 
contribution, which varies as uj/q for small uj/vFq- At 
nonzero temperature, Im^ tp oc T 2 juj. The divergence at 
low oj is not physical, and is a consequence of our having 
neglected the effects of real processes on the matrix ele- 
ment to two-pair states. When such effects are included, 
the divergence will be cut off. This is an example of the 
Landau-Pomeranchuk-Migdal effect . As pointed out 
in Ref. [l(j in the context of neutrino scattering, these 
will become important when uj becomes comparable to 
the quasiparticle lifetime (Eq. 12411 ) which, by manipu- 
lations similar to those used for Imx tp may be written in 
the form 



= 2tt 



(oj + £l ) 2 + (nT) 2 



E 



m*(m*) 2 



(27r) e 



spins 



1*4 



A 2j 



A 4j 



Ar 



o 



K32) 



where the prime on the X^pins indicates that the sum 
over spin for quasiparticle 1 is excluded. In Fig. [21 we 
show a schematic plot of the behavior of Imx f° r small q. 



The quasiparticle scatterin g ampli tude has the follow- 
ing form for neutron matter ( 
on the Fermi surface) 



g ampli 

mm 



) (for quasiparticles 



•A-trtr' (Q: Q j P) -^-scalar ~\~ -^spin cxch.^" ' & ~\~ 
pin— orbit 

i((T + a') ■ q X q' + Aensor^q) + 
^cxch. tensor S(q') + A cm tensor 5(P) + 

A&m. vcctor«(cr - cr') • q x P + 

vector (<T X a') ■ (q' X P), 



(33) 
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FIG. 3: A schematic picture of the behavior of the response 
for small q. The bump at uj < VFq is the contribution from 
single-pair states, the spike is due to the collective mode. The 
straight line represents the behavior of the two-pair response 
for zero temperature, while the dashed line represents the two- 
pair response for nonzero temperature, neglecting the multi- 
ple scattering effect, which is included in the dotted line. 

where |17] 

S(q) = 3cr ■ qer' • q <x • <x', (34) 

and 5(q') and S(P) have the same structure, with q 
replaced by q' respective P. where q = Pi — P3, q' = 
Pi — P4 and P = p! + p 2 = P3 + P4. The amplitudes A 
are functions of the scalar invariants that can be made 
from the vectors q, q' and P. For particles on the Fermi 
surface q, q', P are orthogonal to each other, i. e. qq' = 
0, q-P = and q'P = 0. For asymmetric nuclear matter, 
the Fermi surfaces are different, and in the case when 1 
and 3 are of species i and 2 and 4 are of species j, this 
means that q is orthogonal to both q' and P, but q' is 
not orthogonal to P, which leads to the following extra 
terms in the quasiparticle scattering amplitude [l5j: one 
similar to the spin-orbit interaction: 

so<3' ' + <*') ■ q x P, (35) 
one similar to the difference vector one 

•tfp d iff q' • p»(* - t') • q x p (36) 

and one similar to the quadratic spin-orbit 

^Pquad, so q'-p(<r-qV.p+ 

er ■ Per' ■ q' - jj q' • Per • er') . (37) 

leading to the following scattering amplitude: 
AZ, (q, q', P) = A% aUr + A% in cxch .cr • a> + 
■^pin-orbit^ + t') ' 4 x q' + ^ cnsor 5(q) + 
"^exch. tonsor^'(q ) + 4m tensor'S'(P) + 



^diff. vector*^ - *') ■ 4 X P + 
Across vector^ X a') ■ (q' X P) + 
4'P soq' ' + cr') • q X P + 

•4'p di ff .q' • p ^ - t') • q x p + 
A% quad r . so q' • p • qV • p+ 

cr • Per' ■ q' - | q' • Per • a') , 

(38) 

where the last three terms are zero for i = j. 

We now comment on which terms in the interaction 
are relevant in specific cases. For the density response, 
the total nucleon number density n p + n n , satisfies a lo- 
cal conservation law, from which it follows from Eq. I|3U|) 
that lim^^o linig->o Imx tp = 0. Similar arguments ap- 
ply for the isospin response. For the total spin response, 
lim^^o linig^o Imx p is nonzero, and all terms in the in- 
teraction apart from the scalar and spin-exchange terms 
contribute. For the spin-isospin response, all terms ex- 
cept the scalar term contribute: for this case the spin- 
exchange term contributes because, while it conserves 
spin, it does not conserve spin-isospin. 

VI. CONCLUDING REMARKS 

Landau's theory of normal Fermi liquids is usually re- 
garded as containing only the effects of single-pair states. 
In this paper we have shown how the leading contribu- 
tions to the two-pair response at frequencies low com- 
pared with the Fermi energy and wave numbers small 
compared with the Fermi wave number may be evalu- 
ated in terms of the collision rate for two-quasiparticle 
scattering. The treatment brings out clearly the impor- 
tant role played by conservation laws. Our main result is 
Eq. (|30|l . Our calculation brings out that there are two 
time scales in the response: one is the energy of a single 
pair, vpq and the other is the quasiparticle lifetime. The 
structure on the energy scale v-pq is due to screening and 
to the presence of the intermediate quasiparticle close to 
the energy shell. 

An important result of our analysis is that, for an op- 
erator which does not represent a conserved quantity, the 
matrix element to a two-pair state with excitation energy 
u) varies as 1/lu. Examples of this behavior in specific 
applications have previously been found in work on pair 
bremsstrahlung Refs. 0, 0, 0, 0, |2(j and neutrino 
scattering |l0j . At nonzero temperature the growth of the 
matrix element for uj — > will be cut off by the Landau- 
Pomeranchuk-Migdal (LPM) effect. Neglecting the LPM 
effect, the multipair contribution to Im% at zero tem- 
perature varies as uj and atnonzero temperatures small 
compared with the Fermi energy by T 2 /uj. 

Since completing the above work, we became ac- 
quainted with Ref. 23j in which two-pair contributions 
to response functions were calculated. The calculations 
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were performed for a central interaction, and therefore 
from the general arguments that we have made, at long 
wavelengths one would expect the contribution of mul- 
tipair states to vanish identically. In the calculations of 
Ref. |23| the wavenumber of the response was nonzero, 
and the authors demonstrated that there is a large can- 
cellation between self-energy corrections to the particle 
propagators and vertex corrections: this is due to the 
fact that their interaction conserved both particle num- 
ber and total spin. The contributions to the response 
investigated in the present paper are the ones that dom- 
inate at long wavelengths, and they have their origin is 
completely different: they are du to the noncentral na- 
ture of the nucleon-nucleon interaction. 

For applications to astrophysical situations, it is im- 
portant to extend the long-wavelength results to finite 
wavelengths and higher frequencies. Of particular impor- 
tance for these applications is the evaluation of the two- 
quasiparticle scattering amplitudes. These have recently 
been calculated by Schwenk and Friman [l^j including 
all contributions to second order in the low-momentum 
effective interaction Vi ow fe [2lL 122^ , and it is desirable to 
include higher-order effects. In the case of neutrino pair 
bremsstrahlung, lo > eg, and therefore uj/vpq > c/vp 
which is significantly larger than unity. Consequently, it 



is a good approximation to take the limit vpqjto — 0, as 
has generally been done 0, 0, . In the case of neu- 
trino scattering, effects on the energy scale vpq could be 
important, since the frequencies of interest are less than 
cq. 

Another direction for future work is to extend the cal- 
culations to operators that change the isospin. Examples 
are charged-current weak interactions, which enter in the 
rate of the modified Urea process. 
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APPENDIX: VERTEX CORRECTIONS 

Here we calculate the vertex correction factors A de- 
fined in Eq. Q. In terms of Landau theory, the vertex 
correction is given by 



_ 5e r 
p = 6W 



k p' 



(A.l) 



neutrons and for protons over momenta, we obtain the 
relations 

and 

/ F pp \ frn p f pn ()n n 

arfng = qv p (l + £±-\!^ + , v p A p (0)^-^l(A.8) 



In the limit of no collisions between quasiparticles, the 
linearized quantum kinetic equation reads 0, P-21]: 



_ ( £ p+q/2 _ £ p-q/2)] 



Sni 



(n 



p+q/2 



p-q 



/2)K = o 



where the deviation of the quasiparticle energy from equi- 
librium is given by 



3 P' 



(A.3) 



Here U l is the external potential acting on species i, and 
fp p , is the Landau quasiparticle interaction. Since we 
consider only systems which are not spin polarized in 
equilibrium, we shall, for simplicity, suppress spin in- 
dices. The calculation of the spin-density response is 
completely equivalent to the calculation of the density 
response, apart from the replacement of spin-symmetric 
Landau parameters by the spin-antisymmetric ones. For 
definiteness, let us consider an external potential U n 
which couples to the neutron number. Inserting Eq. 
(|A.3jl into Eq. (|A.2j) gives 



( e F>+q 



/2 



t p-q/2 



)]<K 



in 



p+q/2 



l p-q/2 



)(^+EEOK.') = °- ( A - 4 ) 



These equations are an expression of the conservation 
laws for neutron number and proton number. From Eqs. 
(lA~7)l and (JO), we find 



(A.2) M = 
3 



1 



n uj N»(0)fF x P 



1 + F pp /3 + Ff'73 + [F pp F? n - (F pn ) 2 ]/9 



(A.9) 



and 
Sn p 



1 



w A p (o)/r 



1 + F pp /3 + Ff n / 3 + [^i PP ^i nn - (^D 2 ]/9 

(A.10) 



where 



Rewriting Eq. (|A.4J) as 



H=4(- 5C/i - EE /pp'H')- 
j p' 



(All) 



(A.12) 



We expand 8riL in Legendre polynomials and write 



5nl 



1 n 



p+q/2 



'p-q/2 



N' l (0) e 



p+q/2 p-q/2 



^^(p-q), (A.5) 



where N l (Q) — m*pp/TT 2 is the density of states at the 
Fermi surface. The normalization factor is chosen so that 
the change in the total density of species i is Suq. The 
Landau quasiparticle interaction is also expanded in Leg- 
endre polynomials 



f 

■I U 



E/^(p-p')- 



(A.6) 



with 



(" 



p+q/2 "p-q/2^ 



p+q/2 £ p-q/2 



)) 



(A.13) 



expanding 5n p according to Eq. I)A.5|) and summing over 
momenta, we find two more relations for Suq, Sn p , Snf 



and 8n p : 



K = Xl{-5U» - [F™5nl + -5n\] - 

o qVp 

N»(0)[f p °6n p + ^^-5n p ])(A.U) 
o qVp 



We now solve Eq. (|A.2|I . Since tensor contributions to the 
quasiparticle interaction are generally small compared 
with the scalar and spin exchange terms 0, , we shall 
neglect the tensor terms. In addition, we shall take into 
account only the I = and I = 1 terms in the scalar 
and spin-exchange terms. By summing Eq. I|A.2|) for 



and 



K = X P (-[F pp 6n p + ^--^5n p }- 

r pn 

^(Q)[/rK + 4--p«nf]), (A.15) 
o qVp 
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where Xo 1S 



For a one component system one finds 



Xo 



p 

N l (Q) 



1 dfi (i 

_ x 2 fj, - (uj/v l F q) 



(A.16) 



with fx = p • q. Inserting Eqs. I|A.14|I and (|A.15(I into 
Eqs. ljA.91) and (|A.10(I gives for the response functions 



^ = -xS(i + « p xg)x 



1 + aP X p + a n Xo + « p xga n Xo + Xo^^Xo 

(A.17) 



and 
<5C/ n 



1 + aP X g + a n X§ + a p xga n x8 + xS/3 n /? p Xo ' 

(A.18) 



where 



a p = F pp 



F PP + [ F PP f nn _ (fP")2]/ 3 



q«P y 1 + FfP/3 + Ff n /3 + [F PP F^ - (F pn ) 2 }/9 

(A.19) 



and 

f3 p = N p (0) [/ pn + 



g y t£«$ 1 + FfP/3 + F^/3 + [F^Ff 11 - (FP n ) 2 ]/9_ ' 

(A.20) 

The corresponding results for a n and /3 n are obtained by 
interchanging p and n in Eqs. I|A.19|> and (|A.20|) . If we 
instead use an external potential SU P that couples to the 
proton number we get analogously, 

5n n 

Xg/3 n xg 

1 + aP X ^ + a n Xo + a p Xo« n xS + XoP n P p Xo 

(A.21) 

and 

5n p 



5UP 



-xg(i + « n xg) 

1 + aP X P + a n xS + u p X P a n Xo + XoW p x£ ' 

(A.22) 



su 



-XoiV(O) 



1 + F oXo + {Lo/v F qYF lXo l{l + Fx/3) 



and 



Sni 



upgl + Fi/3 



Sn , 



(A.23) 



(A.24) 



We can now find an expression for the vertex corrections 
from Eq. (|A.1(1 . Expanding Sn p ' a nd f^ , in Legendre 
polynomials, (Eq. (|A.5fl and Eq. l|A.6|l ). and keeping 
only the / = and I = 1 terms gives 



E -, tint, 
eik P 
Jpp' " 



5U3 



3 8W 



(A.25) 



and therefore 



5U« 



= 1 



fS 



np 



(p ' q)w/{qvl)f* } 



1 + FfP/3 + Ff n /3 + [FfFf 
1 



(FD 2 ]/9 



iTinn 
^0 



[ F nn F PP + (F pn ) 2 }/3) 



1 + FfP/3 + Ff n /3 + [F pp F? n - (FP n ) 2 ]/9 



5U n 



5n% 
5U» 



(A.26) 



and 

K 

6U n 



./o 



(p • q) U /(qv%)f? 



l + FfP/3 
1 



Ff n /3 + [F pp F? n - {F pn ) 2 }/9 



ppp 



7VP(0) 
(p-<i)u/(qv»)(FF 



[F „n F PP + (i?P n )2]/ 3 ) 



1 + if p/3 + Ff n /3 + [F pp F[ m - (FP n ) 2 ]/9 



For a one component system this reduces to 



SU n 



8U n ' 
(A.27) 



5U 

which gives 
5e 



Se p 5n .„ _ fx <5ni 

— = 1 + /o 7L7 + (p ' q) y«7' 



(A.28) 



c p 

SU 



= 1- 



p q 



I 



Fxo 



npgl + Fi/3 



1 + [F + (u/vpqfFi/il + Ft/3)} X o ' 

(A.29) 

For uj 3> vf<7, the vertex functions reduce to dei/5Uj = 
Sij. The response to a field that couples to the nucleon 
spin is given by exactly the same expressions, but with 
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the spin-independent Landau parameters Fi replaced by culated here correspond only to the quasiparticle (Lan- 
the spin-dependent ones G/. We note that in the case of dau) contribution, and do not include the part due to 
the spin-dependent response, the response functions cal- multipair intermediate states. 



